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Abstract 

The Fluctuation Relation for a stationary state, kept at constant energy by a deterministic 
thermostat - the Gallavotti-Cohen Theorem — relies on the ergodic properties of the system 
considered. We show that when perturbed by an energy- conserving random noise, the rela- 
tion follows trivially for any system at finite noise amplitude. The time needed to achieve 
stationarity may stay finite as the noise tends to zero, or it may diverge. In the former case 
the Gallavotti-Cohen result is recovered, while in the latter case, the crossover time may 
be computed from the action of 'instanton' orbits that bridge attractors and repellors. We 
suggest that the 'Chaotic Hypothesis' of Gallavotti can thus be reformulated as a matter of 
stochastic stability of the measure in trajectory space. In this form this hypothesis may be 
directly tested. 
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I. INTRODUCTION 



The Fluctuation 



<t T =i£dta(t) 



'. 'heorems are relations satisfied by the distribution of entropy production 



P(-a T ) 



or 



( e - Xr ^) = (e 



-(1-A)tct t 



) VA 



(LI) 



They hold in a series of different settings. The proofs are very simple and general when they 
apply to systems starting from an equilibrium configuration 2 L or in contact with a thermal 
bath involving stochastic noise that guarantees ergodicity 5|, |6) . On the contrary, when the 
system is treated in the stationary state achieved thanks to a deterministic thermostat, the 
derivation is not simple and requires a knowledge of the ergodic properties that is very rarely 
accessible analytically [3|, LZI] . 

In those cases in which the nature of the thermostat is thought to be physically irrelevant, 
it is quite natural to prefer the stochastic baths because they invoke the least hypotheses, 
making the analysis conceptually simpler. On the other hand, the very fact that for determin- 
istic systems the Fluctuation Relation is not automatically satisfied makes it an interesting 
tool to study questions of ergodic theory that are still largely open. 

For a purely deterministic forced system converges to an 'attractor' set, while 

the time-reversed dynamics converges to a 'repellor' set. Loosely speaking, the Gallavotti- 
Cohen theorem is then based on two types of hypotheses: i) that the attractor is sufficiently 
chaotic, and ii) that the attractor and repellor are interwoven fractals with overlapping 
closures. As the forcing is increased, attractor and repellor tend to separate, and even if 

n 

the attractor remains chaotic the Fluctuation Relation no longer holds in general [8J. For 
macroscopic systems at reasonable levels of forcing, one can assume that the hypotheses 
above are effectively (if not always strictly) satisfied, leading to the Chaotic Hypothesis 9]. 
In this paper we recast the content of the previous paragraph in the language of stochastic 



stability 



101 ]: we first show that adding a small, energy-conserving noise the Fluctuation 



Relation is trivially satisfied by any model, and we then analise the limit of noise going 
to zero. Indeed, all the subtelties of the Gallavotti-Cohen theorem and of the Chaotic 
Hypothesis are a result of this limit. The small noise limit has moreover the advantage that 
an approximation scheme (analogous to the semiclassical treatment) becomes exact, so that 
one can analise it in detail. 
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II. THERMOSTATTED DYNAMICS 



We shall consider Hamilton's equations with forcing, a thermostatting mechanism, and 
energy conserving noise 111 ]: 



Qi = Pi 
Pi = - 



dV(q) e , . , , dV(q) 

+ QijVj ~ /i(q) + l{t)Pi = a + QijiVj ~ fj) 



% ' * lJ U Jl ^' ' ,v ^ dq. 
here and in what follows we assume summation convention. 



(II.2) 



rjj are white, independent noises of variance e. 



9ij = ~~ ^ s the projector onto the space tangential to the energy surface. 
/i(q) is a forcing term that depends on the coordinates only. 

l{t) = Multiplying the first of ( 111.21) by g , the second by pi and adding, one 
concludes that energy is conserved. 

The product g%jT]j is rather ill defined because both terms depend on time and rjj is 
discontinuous. The ambiguity is raised by stipulating that it has to be interpreted in 
the Stratonovitch convention 12|, the meaning of which will be made clear below. 



The probability evolves through 121 ] : 



P = -HP 



(II.3) 



where H is the operator: 



d 8V(q) d d r . d . d 
dqi dqi dpt dpi dp { dp. 



■9ij9u 



d_ 

dpi 



(II.4) 



The precise factor ordering in the last term is exactly what we meant by Eq. ( 111.21) being 
in the 'Stratonovitch convention'. In the absence of driving f = it is easy to check that 
H anihilates any function that depends on the phase-space coordinates only through the 

2 

energy E = ^- + V. Hence, the noise respects the microcanonical measure, in that case. 

We shall consider as usual the 'time-reversal transformation' H = [QHQ^ 1 ]* 1 , where 
Q is the operator that reverses momenta: QpiQ^ 1 = —pi and Q^-Q 1 — Under 
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momentum reversal all terms change signs, except the last. Under Hermitean conjugation 
the factor order is reversed, and derivatives change sign. We get: 

H = [QHQ- 1 ]^ = H + a 

» s b»C "l-M = -(«-lb (H.5) 

with 2d is the dimension of phase-space, which can be checked by straightforward calculation. 
We hence have: 

[QHQ- l ] ] = H + a (11.6) 
[Q(H + X^Q- 1 ] ] = H + (1 - X)a (11.7) 

because the entropy production rate a changes sign under time-reversal and conjugation. 

From here it is standard to show that for every finite value of e the fluctuation theo- 
rem holds For completeness, let us outline the proof: given any initial and final 
configurations Xi = (qi,Pi) and Xf = (q/,Pf) the expectation of the total work is 

( e - x $> dt ) = (x f \e-< H+x ^\x t ) 

= ( Xi \e-< H + X ^\x f ) = (xilQ-h-^+^-^Qlxf) (11.8) 

Because of the noise, the evolution is ergodic on the energy surface, and thus the real part of 
the spectrum of H restricted to functions on the energy surface has a gap. Evaluating fill. 81) 
on a basis of eigenf unctions on the energy surface, the 'lowest' right and left eigenf unctions 
(the ones with the eigenvector /x (A) having the lowest real part) (L \ and \R Q ) dominate: 

(x f \e- T ^ H+XT ^\xi) ~ (xf\Ro)(L \xi) e- r ^ (A) 

= (x f \L }(R \ Xl } e-^d-A) (HQ) 

and using the fact that H + Act and H + (1 — X)a have the same spectrum, we conclude that: 

lim l-\n(e- x f° (7dt ) - - hx^ 1 -***? odt ) \ = (11.10) 
t— >oo yj r J 

which is the Laplace transformed version of the Fluctuation Theorem (II. ip . Clearly, the 

whole argument breaks down at strictly zero noise, when the eigenfunctions on the energy 

surface are no longer smooth, and the overlaps (L Q \xf) and (R \xi) may vanish. The problems 

thus may (and will) arise because we have proven the theorem for 

lim lim -j-hipf/ dt a J (11.11) 
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and the Gallavotti-Cohen theorem is for: 



lim lim - InP ( / dt a ) (11.12) 



t— >oo e— *0 T 

In what follows, we shall analise the small-noise limit, starting from a simple example. 

III. A SIMPLE EXAMPLE 

Consider a particle on a two-dimensional space with toroidal boundary conditions, a 
constant field E along the direction 1, and an isokinetic thermostat. The example was 
discussed in Ref. {tJ. We write the velocity vector (p x ,p y ) as: 

p x =pcos9 ; Py = psin9 (III. 13) 

There are two stationary situations: parallel velocity (p x ,Py) = (p, 0), 8 = (stable: the at- 
tractor) and antiparallel velocity (p x ,p y ) = (—p, 0), 9 = it (unstable: the repellor). Equation 
(111.21 ) reads, in this case: 

Epl P X T] X + PyTJy 

Px = E ^ + Vx~Px — 2 

EpxPy PxVx+PyVy frniA\ 

Pv = -2— + r ly-Py -^r^ ( IIL14 ) 

In terms of 8 both equations collapse into a single one: 

9 = -E sin6 + r]g = -JL[-E cos 6] + r] 9 (111.15) 

da 

where we have defined the angular, isotropic noise: rjg = sin Qr\ x — cosdr] y , a white noise 
whose variance is e. The entropy production rate is a = Ecos6. 



Weak noise limit 



Equation ( 1III.15I) can be interpreted as a Langevin process in an effective potential V 
-EcosO, see Fig. [TJ The probability evolves according to: 



d 2 



P 



-H FP P 



(111.16) 



If we are interested in calculating the average e x f adt in the large-time limit we need to 
compute the lowest eigenvalue of 



H x = H FP + XE cos 9 



;iII.17) 
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which can be taken as usual to the Hermitian form 



ei?sin( 



-eEsinO 







E 



e 7^ + T sin 2 # + e£(A--)cos, 



;m.i8) 



and takes the form of a Shroedinger operator in an effective potential V e ff — -f sin 2 9 + e(A — 
\)Ecos9 depicted in the figure Q] The average work per unit time at long times is given 
by the lowest eigenvalue of (lIII.18j) . The Gallavotti-Cohen symmetry is manifest making 
(A - 1/2) -> -(A - 1/2) and 9 -> tt - 0. 




FIG. 1: Top: potential in terms of angle. Bottom: effective potential associated with the 
Shroedinger problem (|III.18|) . A is the subdominant bias proportional to (A — 1/2) that favours 
attractor (A > 0), or repellor (A < 0). 

Now, in the limit of small e the situation is as follows: to leading order in e the ground 
state is doubly degenerate, the eigenfunctions being concentrated in 9 = and 9 = tt, 
respectively - both attractor and repellor are selected by the potential V e ff at this order 
(e -1 ). It is the next-to-leading order A = (2A — l)i?cos# in V e ff (see Fig. [Tj that lifts 
the degeneracy by a quantity of 0(1) in e, thus selecting attractor or repellor, depending 
on the sign of (A — 1/2). Clearly, if (A — 1/2) is of order one, for small e the ground state 
eigenfunction is peaked in only one well, which one depends on the sign of (A — 1/2), and : 

( e ~ A ^> „ e r(|A-l/2|-V2)E (m i9 ) 

where the last term 1/2 comes from the zero-point energy of the minimum. 
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B. Attractor, repellor and time- reversal 



Let us try now a more direct approach (the techiques can be found in 13] ) . We can write 
an expression for the probability of a trajectory in configuration space: 

P[traj) = e~-*S dt [l(*+^) 2 -^"H (m 2Q) 

and we wish to calculate the average 

(e~ Xa ) = [ P(traj) e~ x ^ dt cosd = [ e -\ S dt[\(e+Es\nef+{\-\)Eec OS e\ (m.21) 

J traj J traj 

The trajectory weight in 21 [) is composed of two terms, the leading term corresponding 
to the Lagrangian: 

L — -(9 + Esin9) 2 = -9 2 + -E 2 sin 2 9 — ± E ^1 (m.22) 
4 V ' 4 4 2 dt V ' 

and a subleading term proportional to (A— 1/2). The dominant terms must satisfy Lagrange's 
equations derived from flIII.22j) : these correspond to a problem with inertia in the inverted 
potential —V e ff- The solutions are: 

9 = ±^A£ + E 2 sin 2 9 (111.23) 



where £ is a constant playing the role of an energy associated with the Lagrangian flIII.22j) . 



It is easy to see that the integral in the exponent ( 1III.20I) (the action) can only be finite for 
large times if the trajectory spends most of the time in attractor or repellor, and this is only 
possible if £ = 0. We thus have that solutions are either solutions of the original noiseless 
problem: 

9 = —Esin9 = [-E cos 9] (111.24) 
d9 

which converge to the attractor and have zero action, or 'rare': 

9 = Esin9 = -j^[—E cos 9) (111.2b) 
d9 

taking from attractor to repellor, and these have a finite action equal to E J d9 9 sin 9 = 2E. 
The general trajectory selected by the first term of ( 1111.211) is composed of sejours in the 
attractor and in the repellor, with rapid descents from repellor to attractor (all three with 
zero action), plus the 'instanton' trajectories with action 2E, see Fig. [2l 



8 



repeller 




attraclor 



FIG. 2: An 'instanton' - 'antiinstanton' gas configuration. The average time spent in the attractor 
and repellor determines the entropy production. 

An important point to note is that the trajectories associated with the Lagrangian (1III.22I) 
in the noiseless limit are, as we have mentioned, many more than the trajectories of the orig- 
inal dynamic system in the absence of noise - they include all the non-zero action solutions. 
Also very important is the fact that the dynamics associated with (IIII.22P differs from the 
original dynamics in that the original stable point 9 = is now unstable, due the reversal 
of the potential. These questions will be relevant in the following sections. 

The exact proportion of time spent in attractor r attractor and repellor r repe u or will be 
determined by the second, subdominant term. This term then gives an action ~ eE(X — 
l/2)(T repe u or — T attr actor) ■ The dominant contribution for the A-dependent average is given 
for large times r by the trajectories that maximise 

- \nP{traj) ~ -e(A - l/2){r repeUoT - r attractor ) + N instanton A/ e (III. 26) 

where N instanton is the number of climbings from attractor to repellor. If A = the system is 
unbiased and it will prefer to stay in the attractor. If, on the contrary, we wish to compute 
the large-deviation function with A ^ 0, then the trajectories that dominate will share time 
in attractor and repellor, so that: 

{^repellor ^~attr actor) @t 7 repellor® r repellor T attr actor & attractor (j~attractor ^"repellor)^- 1 (III. 27) 

For the probability of work per unit time a , at small noise level, we get a linear profile as 
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in figure [3] 
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FIG. 3: lnP(a) vs. a for the problem of section II. The profile is linear, and can be seen as the 
coexistence curve of two phases with free-energy ±E. 

which is in agreement with flIII.19p . 

For infinite times the last term in flIII.26j) is negligible, but for times such that instanton 
configurations weigh comparably to the terms proportional to time: A/e ~ In r and we shall 
not see reversals. There are then two possibilities: we may consider a set of initial conditions 
corresponding to an equilibrium configuration, and get the correct proportion of sejour in 
attractor and repellor without needing to perform any activated jumps from attractor to 
repellor. This is the transient fluctuation theorem, which holds for zero noise. On the other 
hand, if we consider any other initial configuration, we need jumps from attractors to repellor 
and back in order to obtain the right proportion of sejours in the large time limit. Clearly, 
this second 'stationary' fluctuation theorem needs lnr ^> A/e. 

Let us note the close analogy to a one-dimensional Ising model E = — J^ 8 SjSj+i — 
h^Si with (3 J — > exp(A/e) the wall energy (where A is the instanton action), and (3h — > 
e(A — 1/2). Although in our case 'climbing' walls cost, and 'descending' ones do not, the 
difference is unimportant because the number of instantons and of antiinstantons differ by at 
most one. Here e enters in the length needed to accomodate many walls, which is necessary in 
the thermodynamic (large- length) limit. The alternative is to consider a properly weighted 
set of boundary conditions, so that no walls are necessary: this is analogous to the transient 
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fluctuation theorem. 



IV. SMALL NOISE LIMIT, GENERAL 

We shall now show that the structure of the previous section is quite general, although 
there are many instanton-antiinstanton trajectories. Consider the evolution equations ( 111.31) 
and (111.41) . We write this as a path integral 

P(q, p) = f e"^ = f e~^ dt (W) (IV.28) 

J [paths] J [paths] 

where the paths have the appropriate initial condition in phase-space, and all possible final 
configurations. 

L = ft (ft - Pi) + piipi + Vi + gijfj) + gijPjgupi 

= Mi+PiPi-H (IV. 29) 

with 

H = qipi - pigijfj - piVi - gijpjgupi (IV.30) 

Here and in what follows we use the notation V: — t^-. The second, subdominant term in 
the exponent comes from the symmetrisation (cfr. the second of (1H.5|) ): 



d 1 

= 2 



\ (IV.31) 



where again a = —id — 1)7. 

We wish to calculate the large-deviation function 

( e X So dt °)=( e -7/<ft[£(q,P)+e(A-§)a] (IV.32) 
J [paths] 

Just as in the previous section, in the limit of small noise we have to look for saddle-point 
solutions only taking into account the first term in the exponent. The equations of motion 
then read: 

; 8H A df r A OVi 

oqi dqi dqi 

Pi = -^ = -ft-M + +M^ 
dH 

ft = ~^r=Pi 
aqi 

OH 

Pi = = -gM + 2pj) - Vi (IV.33) 
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Eqs. ( 1IV.33|) are Hamilton's equations, and they involve twice as many degrees of freedom 
as the original system. In fact, the original phase space variables (qi,Pi) have now become 
the coordinates, and their conjugate momenta are the new variables (qi,Pi). Multiplying the 
third equation by Vi, the fourth by pi and adding, we conclude that all the solutions in the 
extended phase-space of (1IV.33I) are still thermostatted to J2p1 + V = constant. 

Let us now rewrite the action in Lagrangian form. We first put pigijfj = PWugijfj 
where we have used gugij = gij. Next, we separate tangential and normal parts of Pipi = 
PWuPi + P1P1 PiPi/p 2 - We can complete squares and collect the terms, to get the Lagrangian: 

L = ft (ft -Pi)~\ \Pi + Oijfj + V if + 



2 (Pi + 9ijfj + Vi) + gtjpj 



2 

+ VJ Y ipiPi+PiV) (IV.34) 
p A 



Using the last two equations of motion (1IV.33I) . we see that at the saddle point level the last 
two terms cancel. We hence conclude that the action to leading order equals 

A = J dt L(q, q, q) = J dt [q\ + g^ft + Vj] 2 (IV.35) 

where we have substituted pi = <ji everywhere in L. Note that L has an unusual second time 
derivative q\. All relevant trajectories are the saddle points of the action flIV.351) : 

d^dL _ d_di_ di_ _ o V36 

dt 2 dq\ dt d<ji dqt 

The solutions are determined by initial and final values of both q\ and <ji. In particular, 
trajectories of the original system without noise are the ones having zero action: 

q, ■ uul, • (IV.37) 

Furthermore, in order that the action be finite, we must have that q\ + g%jfj + Vi — > for 
large t. As an example, in the next section we shall construct the solutions of Eq. (1IV.36I) 
for the Lorentz gas. 

In a completely chaotic system, we can construct an expansion in terms of (isolated) 



trajectories 
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151 ]. and just as before these will commute between attractor(s) and repel- 
lor(s) - defined as dynamically stable and unstable solutions of the equations of motion of 
the original system - where they will spend essentially all the time, since otherwise their 
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action would be infinite. Because the only difference between attractors and repellors is, by 
definition, their stability, the next-to-leading contribution in the action 



A e = [ A - ^\ I dt a + \ [ dt dt' 



2J J 2 J 5 qi (t)5 qj (t') 



6qi(t)S qj (t') (IV.38) 

traj 



will decide the length of the sejours in each for each A. 

In this setting, a necessary condition for the Gallavotti-Cohen theorem (i.e. the Fluctua- 
tion relation at zero noise) to hold, is that there be orbits with zero action bridging attractor 
and repellor. 

A. Variational solutions: the zero-action situation. 

Given a point on the attractor (q a ,p a = q a ) and a point on the repellor (q r ,p r = q r ) 
and two times ti and tf we can obtain the solution going from (q a ,q a ) at tj to {q r ,q r ) at 
tj as follows. Starting from any curve (q(t),q(t)), contained within the energy surface and 
having the correct endpoints, its action will be an upper bound for the true, minimal action. 
Deforming the curve so as to minimize the action (1IV.35|) we reach a solution. The arrival 
time may become a variational parameter itself. 

If the lowest action A m j n associated with any trajectory joining any two points in the 
attractor and repellor is finite, we know that the system will necessarily need a time at least 
t ~ e A m i„/e ex ] ore them. This can be simply seen by looking at the sum over paths 
(1IV.32P as a partition function in a space of trajectories at inverse temperature e, with A m j n 



its 'ground state energy'. 

Clearly, for the Gallavotti-Cohen theorem to apply, a necessary condition is that A m ,j n = 
0. For this to happen, attractor and repellor have to be interwoven so that just advancing 
from (q a , q a ) along a zero-action trajectory satisfying flIV.37j) we reach a point that is as close 



as we wish to some (q r , q r ), providing a zero-action 'bridge'. However, the condition A min = 
only implies that the time for bridging attractor and repellor grows in a subexponential — for 
example power law r ~ — way, and is not necessarily finite as e — > 0. One can imagine 
such a situation arising in a case in which attractor and repellor have closures of dimension 
d\, but these intersect in a manifold of dimension di < d%. Without noise a trajectory might 
have vanishing probability of leaving the attractor and entering the repellor, since the time 
spent in the overlapping region would be negligible; the addition of noise will coarsen both 
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attractor, repellor and intersection by introducing a width, hence making the passage more 
likely. 

V. INSTANTONS FOR THE LORENTZ GAS 



We shall consider here a Lorentz gas [la], a system with a single thermostatted particle 
under the action of an electric field as in section II, but with in addition obstacles where the 
particle bounces, introducing chaos The Lagrangian (1III.22I) reads 

I = ^(9 + Esin9) 2 =p e 9-H 

H = E = 2 Pe (po - E sin 6) (V.39) 
with 8 defined in (1III.23D and 

p e = -(9 + Esm9) (V.40) 

The motion is punctuated by bounces on the wall. Between bounces, the trajectories satisfy 
(1ITL23|) : 



d9 = ±^4£ + E 2 sin 2 9dt 

, , cos 9 cos 9 

ax = cos9dt = — : — da = + - riR 
9 VE 2 sin 2 9 + 4£ 

dy = sin Mi = h ^-d9 = ± - Sm ° =d9 (V.41) 
9 VE 2 sin 2 9 + AS 

which yields the parametric equations for the trajectory segment. One also has that 

pe = —2E cos 9 po = 2E pe x 
Pe(t) = Pe(t )e 2Ex \° (V.42) 

where t a is the initial time of the segment. The action for the trajectory segment reads: 

A = [dt p 2 e (t) = p 2 e (t ) [ dt e 4E < (V.43) 



where x\ l = J* dt x is the total distance along x, without subtracting the windings if the 
space is periodic. It is positive if the particle moved downhill, and negative if it moves 
contrary to the force. 

Next, we need to consider bounces. Because we are solving the instanton dynamics, rather 
than the original one, we could wonder what happens during a bounce. To make this point 
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clear, in the Appendix we compute a bounce by considering the wall as a limit of continuous 
potentials. The result is that: i) the reflection law holds, ii) the contribution of the bounce 
to the action vanishes in the hard wall limit, and in) pg does not change from just before to 
just after the bounce. Instead, £ is not conserved in bounces, since: 

£\^Ze=^E sm9\fXre (V.44) 

which is nonzero in general, except when pg = 0. 

Segments of trajectories belonging to the attractor and repellor of the original noiseless 
dynamics have zero action and pg(t) = £ = 0: this is preserved by both trejectories and 
bounces at all times (cfr. Eqs. flV.43j) and flV.441) ). 

Trajectories have a starting point (x,y,6,6). Thanks to bounces, they move uphill and 
downhill with respect to the field. Those that diffuse uphill have a smaller and smaller value 
of pg ( Eq. (1V.43I) ). Furthermore, as pg — > also £ — ► (although £ can still increase 
somewhat in a bounce). Hence, a trajectory that had bounces allowing it to diffuse a long 
way uphill becomes more and more closely a zero-action trajectory. Moreover, it is (or 
rather, it shadows) a trajectory belonging to the repellor, since trajectories in the attractor 
move on average downhill. 

Suppose next that the starting point is a very small perturbation with pg ^ of a point 
on the attractor. At the beginning the point will continue to move on average downhill. 
However, this implies that the value of pg grows exponentially ( Eq. (IV.43I) ). In other words, 
in the full dynamics with instantons, true attractor trajectories are unstable, a general fact 
already mentioned in the previous sections. After some time, these trajectories will bounce 
uphill and downhill. A few of these perturbations will diffuse a large amount uphill, and 
as described above those will then be very close to the repellor. Hence, we have identified 
the instantons interpolating between attractor and repellor as the (rare) trajectories that 
have initial conditions close to the attractor and such that they keep diffusing on average 
upwards. Their action is just the integral of the exponential of the uphill distance - a finite 
quantity. 

As is typical with instantons, we need some sort of 'shooting' method to find those 
solutions that actually end up in (or shadowing an orbit in) the repellor, otherwise a generic 
solution will just evolve away. Secondly, we see that there are many of these solutions, and 
not essentially one as in the case without obstacles. 
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VI. CONCLUSIONS 

We have recasted the problem of the deterministic Fluctuation Theorem as the vanish- 
ing noise limit of stochastic one. Large-deviation functions in this limit need not coincide 
with the zero-noise result, but one can always postulate that this is so in a practical case, 
thus making an alternative Chaotic Hypothesis |9|j: Considering the energy-conserving 
noise described above (and more generally noise respecting all the constants of motion of 
the dynamics), a reversible many-particle system in a stationary state can be regarded as 
stochastically stable for the purposes of computing probabilities over trajectories of macro- 
scopic observables. The advantage of this formulation is that it is immediately testable in a 
numerical (and perhaps also in a real) experiment. 

We have argued here that if there are orbits with zero action that bridge attractor and 
repellor, then the time needed in order that the Fluctuation Relation is satisfied is subex- 
ponential in the noise elnr — > 0. We have, however, fallen short of proving the Gallavotti- 
Cohen theorem for that case, since we argue that subexponential yet slowly divergent growth 
of r as e — > is likely to arise, perhaps in systems in which the closures of attractor and 
repellor are simultaneously dense in a manifold of low dimension. It seems an interesting 
question to obtain estimates of passage times in terms of the noise (or temperature) in such 
cases. 
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VII. APPENDIX. INSTANTONS AND HARD WALLS 

To model the bounce, we consider walls as regions with a constant repulsive force (per- 
pendicular to the wall), and then take the limit of large force. We shall rotate the axis to 
a new axes (x',y') so that outside the wall region the field is (E\, E 2 ) and inside the wall 
region the field is (E w , 0), E w — > 00 as shown in the figure. 




FIG. 4: Bounce trajectory in a wall region. 



In general, in a region of fields (Ex, E 2 ) the Lagrangian reads 

L = ^(9 + E lS m9 - E 2 cos6) 2 
n = 2p e (p e - E 1 sin 9 + E 2 cos 6) 



(VII.45) 



with 



1 



£!sin# - E 2 cos9) 



(VII.46) 



The forces Ei jump when the particle crosses the wall. Across this jump, 9 is continuous, 
but 9 jumps. To calculate this, we write the equations of motion: 

2^- (9 + E 1 sin 9 - E 2 cos 9) - 2{9 + E 1 sin 9 - E 2 cos 9) {E x cos 9 + E 2 sin 9) = (VII.47) 



which implies that in regions of constant fields: 
^\n(9 + E 1 sin9 - E 2 cos9) = 



In 



9 + E l sin9-E 2 cos9 



2{E 1 x + E 2 y) 

2(E 1 x + E 2 y)\ t = 2x.E| 



(VII.48) 
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On entering the wall region, equation (1VII.47|) implies that 

p e (f ) = - (0 + E x sin B - E 2 cos 6) ( VII.49) 

is continuous across discontinuities of the potential. Instead, the numerical value of H = £ 
jumps when the trajectory bounces, except if pe = 0. 




FIG. 5: The angle evolves in a potential ~ —Esin 2 9, entrance and exit angle a symmetrically 
disposed, since the distance along x moved is the integral of cos 6>(i), which is exactly zero if 
tt/2 - in = 9 out - it/2. 



Within the wall region: 
The conserved quantity is: 
and we have: 



L = -(d-E w sind) 2 (VII.50) 



£ = 0^E^ sin 2 6 (VII.51) 



9 = ±\JEl sin 2 fl + 4£ (VII.52) 

This is shown in figure and it is clear that the reflection law will hold, since the trajectory 
inside the wall region is symmetric. In terms of the coordinates normal and parallel to the 
wall: 

, , _ , cos# cos 9 

ax = cos 9at = — = — a9 = — rlH 

9 ^JEI sin 2 9 + AS 

dy' = sin9dt = -^-d9 = - =d9 (VII.53) 

9 ^/Elsm 2 9 + 4S 
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We have x! = cos 9 and (1VII48I) reads: 

(0 - E w sin 0) t = (6-E w sin 0) t=o e" 2 ^' (VII.54) 

It is then clear that, in the limit of strong E w , the time spent in the wall region tends to 
zero, while the integrand in the action e AEwX ' is itself of order one (since the penetration 
x' ~ 0(1/E W )). We conclude that the action during the bounce tends to zero as \E W \ — > oo. 
Note that the argument is valid for quite general bounces. 
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